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^-j. , ABSTRACT 

■ The identification of two new Planetary Nebulae in the Sagittarius Dwarf Spheroidal Galaxy 
<~) (Sgr) is presented. This brings the total number to four. Both new PNe were previously classi- 
\Q , fied as Galactic objects. The first, StWr 2-2 1 , belongs to the main body of Sgr, from its velocity 

■ and location. The second, the halo PN BoBn 1, has a location, distance and velocity in agree- 
ment with the leading tidal tail of Sgr. We estimate that 10 per cent of the Galactic halo consists 
of Sgr debris. The specific frequency of PNe indicates a total luminosity of Sgr, including its 

I ' tidal tails, of My = -14.1. StWr2-21 shows a high abundance of [O/H] = -0.23, which con- 

firms the high-metallicity population in Sgr uncovered by Bonaficio et al. (2004). The steep 
metallicity-age gradient in Sgr is due to ISM removal during the Galactic plane passages, 
ISM reformation due to stellar mass loss, and possibly accretion of metal-enriched gas from 
our Galaxy. The ISM re-formation rate of Sgr, from stellar mass loss, is 5 x 10~ 4 M Q yr _1 , 
amounting to ~ 10 6 M per orbital period. 

HST images of three of the PNe reveal well-developed bipolar morphologies, and provide 
clear detections of the central stars. All three stars with deep spectra show WR-lines, suggest- 
ing that the progenitor mass and metallicity determines whether a PN central star develops a 
WR spectrum. One Sgr PN belongs to the class of IR-[WC] stars. Expansion velocities are 
determined for three nebulae. Comparison with hydrodynamical models indicates an initial 
density profile of p oc r~ 3 . This is evidence for increasing mass-loss rates on the AGB. Peak 
mass-loss rates are indicated of ~ 10~ 4 M Q yr _1 . 

The IR-[WC] PN, He 2-436, provides the sole direct detection of dust in a dwarf 
spheroidal galaxy, to date. 

Key words: Galaxies: individual: Sagittarius dwarf spheroidal; Planetary nebulae: extra- 
galactic; Stars: mass loss' Stars: Abundances 
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1 INTRODUCTION 

Planetary Nebulae (PNe) in Local Group dwarf galaxies provide 
valuable indicators of the star formation history. They are espe- 
cially important for the age-metallicity relation, since they allow 
us to break the age-metallicity degeneracy inherent in colour- 
magnitude diagrams. Their abundances can provide stellar nucle- 
osynthesis information for metallicities often very different from 
local, Galactic PNe. They are also important for studying PN evo- 
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lution itself. Knowing the age and metallicity of the stellar popu- 
lation also allows one to study the formation and evolution of the 
PNe as function of these properties. In particular, PNe form dur- 
ing the catastrophic mass loss which terminates the Asymptotic 
Giant Branch (AGB) evolution. T he de pendence of the mass loss 
on metallicity is not understood lZiilstrall2004l : Ivan Loonll200d: 
Ivan Loon et aljEo05l) . and PNe in systems with well-determined 
stellar populations are crucial for studying this important problem. 

The Sagittarius dwarf galaxy (Sgr) is our nearest surviving 
neighbour, at a distance of 25 kpc, less than half that of the Large 
Magellanic Cloud. Its location behind the Galactic Bulge con- 
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tributed to its late discovery llbata et aljfliffll) . The present mass 
and luminosity are low, within the range of dwarf spheroidals. The 
abundances range from SMC-like to LMC-like for different star- 
formation epochs (e.g. Lavden & Saraiedini 200G; Alard 2001). 
Sgr is strongly disrupted by its interaction with the Milky Way 
(Putman et al] |2004 iHelmi & White! 1200 lb iDohm-Palmer et alJ 
120011) . and was initially a substantial galaxy, in accordance with its 
relatively high metalicity. Some distant carbon stars in the G alac- 
tic ha lo may derive from the tidal streams of Sgr iMauron et alJ 
2005). This has also been suggested for distant halo globular clus- 
ters iBellazzini et alj2 003). What fraction of the Galactic halo was 
strippped from S gr is not kno wn. 

IZiilstra & Walshl 1 19961) discovered that two previously cat- 
alogued (Galactic) planetary nebulae were located in the Sagit- 
tarius dwarf galaxy on the basis of their location and radial ve- 
locities. They are thus the closest extra-galactic PN. A detailed 
analys is based on gro und-based spectra and radio continuum data 
iDudz iak et alj |2000) showed that the two nebulae have the same 
stellar progenitor mass of 1.2Mq and almost identical progeni- 
tor light element abundances. One of the nebulae, Wray 16-423, 
underwent PN ejection about 1500 yr previous to its twin He 2- 
436. A differential abundance analysis revealed th e first conclusive 
evide nce for third-dredge-up oxygen enrichment IPeauignot et alJ 
2000). That the two PN have such closely matching ages and abun- 
dances of non-dredged-up elements, pointed very strongly to their 
belonging to the same, short-lived star-formation episode. 

In this paper we present a new member of the Sgr PN popula- 
tion, arising from a recent, higher-metallicity star-formation epoch. 
One further Galactic halo PN is also found to be a Sgr member, 
located in the leading tidal arm. This object traces an older, much 
lower metallicity population. 

Section 2 discusses the PNe membership of Sgr. In Section 3 
we present new HST images for three PNe, and deep spectroscopy 
of the new high-metallicity PN. Kinematical models are discussed 
in Section 5. The central stars and the nebular morphologies are dis- 
cussed in Section 6. Abundances and implications for the evolution 
of the Sagittarius dwarf galaxy are discussed in Section 7. 



2 MEMBERSHIP 

2.1 A Sgr-core member: StWr2-21 

The main body of Sgr lies behind the Bulge, at a Galactic latitude of 
minus 15-20 degrees. The distance is between 24 ± 2kpc l Alard 
1 19961) and d = 26.3 ± 1.8 kpc iMonaco et al] |2004). Confusion 
with the Bulge foreground population is important, but there are in 
fact relatively few Bulge PN more than 15 degrees from the Galac- 
tic Centre. The two known members, Wray 16-243 and He 2-436, 
were identified as located in the main body of Sgr based on their 
identical radial velocities IZiilstra & Walsh 1996). Another object 
located in the same region (Hb 8) was shown to have a very dif- 
ferent velocity and was classified as a foreground object. Further 
PNe in Sgr, wrongly classified as Galactic Bulge PN, are a distinct 
possibility: many Bulge PNe have no known velocity, and there 
are suggestions that the Sagittarius galaxy could have had a total 
mass > 10 9 M Q Ijiang & BinnevlEoQcl) . We therefore conducted a 
search through the literature for catalogued PN in the direction of 
the Sagittarius galaxy, and to test for a radial velocity compatible 
with Sgr membership. 

A little studied PN, StWr 2-21 (PNG005.2-18.6), was iden- 
tified situated within 2° of Wray 16-423. Its velocity was mea- 
sured using a newly obtained WHT spectrum (see Section |4~H . 



The mean measured heliocentric radial velocity of the H/3 and 
[O III] 4959,5007A lines was 129±1 kms" 1 (error on mean). This 
agrees to within 5 km s _1 with the other two Sgr PNe, leaving little 
doubt about the association with Sgr. 

All three identified members are located within the southern 
extension of Sgr, in the same area where its globular clusters are 
found. There is no PN candidate in the core region of Sgr cen- 
tred on M5 4 (centre I , b — 5.57, —14.17, core radius 234 arcmin: 
iMaiewski et aljl2003l) . The only catalogued PN in this region, He 
2-418, has an acceptable velocity of + 1 15 km s _1 but a diameter of 
14 arc sec I Ruffl e etall2004 which is too large for a bright PN at 
the distance of Sgr. Assuming that this region has been surveyed to 
uniform depth, the lack of bright PNe towards the assumed bound 
core region of Sgr is noticable. There is no evidence of high extinc- 
tion in this direction. 

IWalsh et alJ < 19971) suggested that PRMG 1 
(PNG006.0-41.9) was a possible Sagittarius PN. Its abun- 
dance is similar to those of the two Sagittarius PNe I Dudzi aTTet alJ 

2000) . We can now rule out this proposed association. Observa- 
tions with the NTT EMMI spectrometer of this 4-arcsec nebula 
showed that its heliocentric velocity is +20±2 kms -1 (LSR veloc- 
ity +24kms~ 1 ). The Sagittarius stream 20 degrees fro m the main 
body is thought to have a velocity close to 100 kms" 1 IQm taejLalJ 

2001) which does not favour membership. iHelmi & WhitdfcOOll) 
predict velocities at this position closer to 50 kms -1 .) Recent 
mappings show that the tidal plane of Sgr does not pass through 
this position iMaiewski et al. 2003J). We conclude that PRMG 1 is 
a foreground object. 

2.2 A possible leading-tail member: BoBn 1 

We also searched along the plane of the tidal tails I M aiewski et alJ 
l2003tlNewberg et all2003l) . This identified BoBn 1, a well-studied 
halo PN, as a cand idate. Its heliocentric velocity of +174 kms -1 
(Wright et al. 2005) is opposite to the trailing tail at this location 
IMaie wski et al.12 004). However, the positional agreement is suffi- 
ciently close (and halo PNe are sufficiently rare) that the question 
may be asked whether BoBn 1 could b e par t of a s eparate tail. 

A dynamical simulation by lLaw et al J 120051) shows the vari- 
ous streams of debris, separated into trailing and leading tails, and 
associated with different tidal dislocations events. The streams are 
always located in the same orbital plane but are moving in differ- 
ent directions and are at different distances. BoBn 1 is located in 
a direction where the leading and trailing streams cross. At this 
location the leading stream is moving outwards, along the line of 
sight. The particular velocity is model dependent, but the velocity 
of BoBn 1 is within the range of the predictions. The heliocentric 
distance to the leading stream at this location is around 20 kp c. The 
distance to B oBn 1 is quoted as 22.5 kpc lHawlev & Miller! 1 97 8h 
and 16.5 kpc iHenrv et all20 04). (Distances to PNe are difficult to 
ascertain, and the quoted distance is uncertain by at least 30 per 
cent.) The agreement in position, velocity and distance makes an 
association of BoBn 1 and the Sgr leading tail a likely one. 

The leadin g tail is observa tionally not as well constrained as 
the trailing tail lLaw et al 120 05). The identification of further stars 
in the Sgr leading tail will be needed to remove the model depen- 
dence in the proposed association.The accurate velocity of BoBn 1 
may also be useful in constraining the Sgr tidal models. 

No other catalogued PNe were identified as potential mem- 
bers. Other newly discovered PNe in the Sgr core region 
lParke ~t alj|2003l) are extended and faint and certainly not Sgr 
members. The list of possible PNe in the ESO-Strasbourg catalogue 
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Table 1. The Sgr PNe. For the first three objects, stellar magnitude, nebular flux, and diameters are derived from the HST images presented here. Data for 
BoBn 1 are from Wright et al. 1 2005) but the position was redetermined. Stellar classification is discussed below. 



PN 



RA DEC 
(J2000) 



(km s 



mv (star) 
_1 ) (mag) 



diameter star 



logF(H a ) 
(ergcm~ 2 s~ 1 ) (arcsec) 



He2-436 4.8 -22.7 19 3206.70 -34 1257.4 +133 17.25 ±0.1 -11.74 

StWr2-21 5.2 -18.6 19 14 23.35 -32 34 16.6 +129 19.7 ±0.1 -12.26 

Wray 16-423 6.8 -19.8 19 22 10.63 -31 30 38.7 +133 18.85 ± 0.2 -11.55 

BoBn 1 108.4 -76.1 0037 16.03 -134258.5 +174 -12.425 
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of Planetary Nebulae (Acker et al. 1992) contains three further ob- 
jects of interest (K 2-4, Y-C 43, ESO 340-34), but all are nebulae of 
large radius. They are possibly misclassified galaxies. 
The four identified Sgr PNe are listed in TableQ 



3 HST IMAGING 

We obtained SNAPshot WFPC2 images in three different filters 
of the three PNe in the main body of Sgr. We used a sequence of 
two cosmic-ray-split Ha exposures, one V-band exposure and two 
[O III] exposures. For He2-436 and Wray 16-423, exposures times 
are 100 seconds for each Ha exposure, 60 sec for the V-band and 
80 seconds for each [Olll] exposure. The faintest object, StWr2- 
21, was observed for 2 times 140 sec in Ha and [O III]. The objects 
were placed on the CCD of the Planetary Camera (PC). The pixel 
scale is 0.0455". Pipe-line reduced images were retrieved. 

The calibration o f the three filters us ed (F502N, F547M, 
F656N) is described in lO'Dell & Doil Jl999l) . The filter response 
is taken from the WFPC2 handbook. The F656N filter is centred on 
Ha wth a peak transmission above 70%. It has a ~ 10% response 
at the [N II] 6548A line and a few per cent at the 6584A[N II] line. 
The He I line at 6560A can also contribute but this line is weak in 
all cases. The F502N filter has a 60% transmission at the [O III] 
line without significant pollution from other lines. The F547M is 
used as the continuum filter, and is best suited to detect the central 
star, and other possible confusing stars-it is little affected by strong 
emission lines. Nebular bound-free continuum dominates the neb- 
ular flux in this filter. 

Photometric calibration was done with the calibration con- 
stants (PHOTFLAM and the Vega-based zeropoint) listed in the 
WFPC2 handbook. The central stars of all three nebulae are de- 
tected in the F547M images. The Ha image was scaled and sub- 
tracted from the F547 image, to remove the nebular continuum. 
The stellar magnitudes, converted to Johnson V, are listed in Table 
□ They were measured in a circular apertures of 7-1 1 pixels radius, 
using aperture corrections shown in the WFPC2 handbook for the 
V-band filter and the PC. The accuracy is 0.1-0.2 mag, limited by 
the faintness of the stars. In addition to the magnitudes, TableQalso 
lists the line fluxes measured from the images. 

Fig. shows the resulting images. Here the images were ro- 
tated to the same orientation, while rebinning to keep the same 
square pixel size. All images are shown to the same spatial scale. 
The images reveal bipolar structures in all three PNe. The large 
difference in size is clear. The peak intensity drops by a factor 
of roughly 100 between He 2-436 and StWr2-21 (in the depicted 
images the fluxes have been scaled to the same peak for the [O 
III] images). He 2-436 is very compact. The inner radius is only 
just resolved, as shown by the two peaks separated by a single 
pixel. It is significantly elongated. Wray 16-423 shows a compli- 
cated morphology, with an inner bipolar structure surrounded by 



a second, fainter structure rotated at an angle relative to the inner 
countours. The outer contours show elongation. Finally, StWr2-21 
is the largest nebula, and shows an elliptical outer structure and a 
bipolar inner structure. There is some evidence for a wind-blown 
bubble towards the East which in [O III] appears as a closed outer 
shell. As for Wray 16-423, a rotation between the inner and outer 
structures is apparent. 

Fig-IUshows radial intensity profiles along the minor axis and 
along the major axis (90 degrees rotated). These were obtained by 
rotating (and rebinning to the original square pixel size) to align the 
major axis with the image axes, and averaging three adjacent rows 
or columns through the geometric centroid of the source. The V- 
band profiles show the central star, although in the case of He2-436 
only marginally. 



4 SPECTROSCOPIC OBSERVATIONS 

4.1 Low resolution spectroscopy of StWr 2-21 

We observed StWr2-21 with the ISIS spectrometer on the 4.2m 
WHT La Palma telescope on 1999 July 19. The holographic 
H2400B grating was used, and the EEV12 CCD detector with a 
pixel scale of 0.1 lA and a coverage of 445 A. A slit width of 
1 .5 arcsec was employed giving a line width of 5 pixels and a result- 
ing spectral resolution of 0.53A. A single exposure of 900s was ob- 
tained with the slit at a parallactic angle; the seeing was measured 
at about 2 arcsec. The spectrum was centred at 4800A allowing the 
He H 4686A, H/3 and the [O III]4959,5007A lines to be studied. A 
Cu-Ar arc lamp exposure was employed for wavelength calibration 
and a 5th-order fit to the line positions was employed (r.m.s. on the 
fit <0.0lA). An exposure of the spectrophotometric standard star 
BD+33°2642 JOkelll990h was made with a broad (5.0 arcsec) slit 
and used to establish the absolute flux calibration for the StWr2-21 
spectrum. 

A deep, low-resolution spectrum with larger wavelengt h cov- 
erage was obtained with EFOSC2 on the ESO 3.6m telescope IPatatl 
1990). An acquisition image was obtained with an [O III] filter 
(centre 5004A, width 56A; ESO Filter #689) and spectra at a paral- 
lactic angle with a 1 .2 arcsec slit were obtained over the wavelength 
ranges 3050-6100 (grism #03) and 4090-7520A (grism #04). The 
seeing was 1.6 arcsec measured from the acquisition image. The 
exposure time was 1800s, split into two sub-exposures; separate 
300s exposures were also obtained to provide all strong lines un- 
saturated. The measured spectral resolution was 1 1 A for the bluer 
spectrum and 12A for the redder spectrum. Broad-slit observations 
of the standard star EG274 lOkdll990l) were obtained with both 
gratings for flux calibration. 

The spectral images were bias-subtracted and flat-fielded by 
exposures to a continuum lamp illuminating the inside of the dome. 
Wavelength calibration was applied using the exposures to a Cu- 
Ar arc lamp fit by a third order polynomial. The spectrum of the 
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Figure 1. Contour plots in [O III], V-band and Ha of the three Sgr PNe. The peak fluxes are scaled to give approximately the same peak value for all [O III] 
images. A small constant was subtracted from the StWr 2-21 images to suppress the lowest contour. Contour levels are at 0.5,1,2,4,8,12,16,24,28,32,36,40 
scaled counts. At the distance of Sgr, 1 arcsec corresponds to 3.9 10 17 cm 



PN was extracted over the extent of the emission, sky subtracted 
and flux calibration applied from the spectra of EG274. Table [2] 
presents the observed line fluxes determined by interactively fitting 
Gaussians to the lines for the longer exposures, with substitution of 
the fluxes for the saturated lines from the shorter exposures ([O III] 
doublet and Ha). Errors on the line fluxes, from propagation of the 
photon counts on the lines and the error in fitting the underlying 
continuum, are also listed in the Table. 

The extinction is zero within the errors, as derived from the hy- 
drogen line ratios. For comparison, the extinction towards Wray 16- 
423 is E(B - V) = 0.14 jWalsh et alll997l) . consistent with the 
forground extinction towards Sgr. He 2-436 has higher extinction 
due to dust within the nebula. 



4.2 Echelle observations 

He 2-436 and StWr 2-21 were observed with the ESO NTT tele- 
scope during June 2001, as part of a la rger survey of PNe towards 
the Galactic Bulge <GesickietaT]E 006). The echelle spectra were 



obtained using the red arm of the combined imager/spectrograph 
EMMI. Echelle grating 14 was used with cross disperser 3. The 
wavelength coverage is approximately 4300-8 lOOA. The slit width 
was 1", providing a resolution of 60000. The slit length is limited 
to 3 arcsec by the need to avoid order overlap. The seeing (worse 
than 1 arcsec) did not allow us to use the spatial resolution inherent 
in the data, and we summed the spectrum over the slit length. 

Wavelength calibration was done using a ThAr lamp. The 
data was flat-fielded and corrected for the response function us- 
ing a standard star. No absolute flux calibration was attempted. 
The exposure time of 120 sec aimed at obtaining good line pro- 
files for the strongest lines in the spectrum only. Line profiles used 
to obtain velocity fields were determined for four lines. For Wray 
16-423, similar but muc h deeper observations were published in 
iGesicki & ZiilstrJ 120031) . The data reduction used here is as de- 
scribed in that paper. 

The line ratios from these data are listed in Table |2| In the 
following sections we describe models based on these new data. 
These are combined with existing models for the previously known 
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Figure 2. Cross sections of the three nebulae observed with HST. The profiles are summed over three adjacent rows. Dashed lines are profiles along the minor 
axis (typically showing the two peaks) and continuous lines are along the major axis. The insets show the inner 0.2 arcsec of He 2-436, showing the flattening 
in [O III] and Ho, compared with the central peak seen in V which is due to the central star. 



Sgr PNe. No new data was obtained for BoBn 1, and this object is 
not included. 



5 EXPANSION VELOCITIES 
5.1 Description of the models 

The internal kinematics of PNe are dominated by expansion. The 
expansion velocity together with the physical size of the nebula 
gives a measure of the elapsed time since the ejection of the neb- 
ula, the dynamical age. This can be compared to the thermal age, 
defined by the temperature of the central star and its rate of temper- 
ature increase. The latter is a sensitive function of the mass of the 
central star. The dynamical age together with the stellar tempera- 
ture can therefore be used to derive the stellar mass. 

The internal kinematics are domin ated by th e expan sion of 
the nebulae. Hydrodyamical models (e.g. Schonberner et al. 2005) 
predict that the expansion velocity is not uniform, but that differ- 
ent layers (at different radii) expand at different rates. The echelle 
spectra clearly show the broadened lines due to the expansion, but 
determination of the internal velocity field is needed to define an 
overal expansion velocit y. The mass-weighted expansion velocity 
iGesicki & Ziilstr J2000h is used for calculating a dynamical age. 

We will fo ll ow t he usual procedure, as outlined in 
IGesicki & Ziilstr J 120031) . Here we first use a photoionization 



model to determine where within the nebula a certain line origi- 
nates, and use the measured line profiles to fit a velocity field as 
function of radius. The model assumes a spherical symmetric neb- 
ula. The models we use (the Torun models) do not currently imple- 
ment asymmetric nebulae, and neither do we have sufficient obser- 
vational constraints to fit such models. (For instance, asymmetric 
velocity models require knowledge of the angle between the polar 
axis and the plane of the sky, which would require spatially re- 
solved echelle spectra. Our spectra lack spatial information.) The 
Torun models also approximate the star with a black body. How- 
ever, the resulting expansion velocity is found to be relatively ro- 
bust to these assumptions. 

The computer model first solves for the photo-ionization equi- 
librium, where we attempt to reproduce the observed line intensi- 
ties and the observed radial surface brightness distribution. (This 
photo-ionization model is by necessity different from that used in 
section [5] and is not suitable for abundance determinations.) We 
compare the models with both major and minor axis cross sections 
of the nebulae, but do not attempt to fit these individually. The line 
profiles are calculated, taking into account the slit aperture (which 
may not cover the full nebula) and the seeing (which scatters light 
rays from outside the aperture into the slit). The velocity field is 
adjusted to fit the line profiles of all lines simultaneously. In the 
present analysis the computer codes are adapted to work with an 
optimization routine. We applied the publicly available "PIKAIA" 
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Figure 3. Model fits to the line profiles of He 2-436. Left: parameters of the ionization model, showing the radial distribution of the including the density, 
electron temperature, velocity field, and the ionic distribution. Middle: radial surface brightness profiles: model (solid lines) versus the observed major and 
minor axis profiles from the HST images, fig. refprofiles. Right: Fitted line profiles. 




Figure 4. Model fits to the line profiles of StWr 2-21. details as in Fig.l3l 
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Table 2. StWr2-21 Emission-line flux ratios as measured from the low- 
dispersion spectra, relative to I(H/3)=100.0.. 



Kesi A 


Identifi c ation 


wds. nux 


error 


uds. nux 


(A) 










3425 


[Ne V] 


1.73 


0.23 




3440 


OIII 


10.16 


0.19 




3705 


HI 


0.59 


0.08 




3726 


[Oil] 


34.36 


0.16 




3753 


OIII 


3.90 


0.07 




3769 


HI 


3.31 


0.09 




3796 


HI 


4.70 


0.07 




3814 


Hell 


1.51 


0.06 




3835 


HI 


6.67 


0.09 




3869 


[Ne ffl] 


70.97 


0.20 




3889 


H I +He I 


16.69 


0.11 




3922 


He II 


0.40 


0.11 




3969 


H I + [Ne III] 


35.89 


0.15 




4026 


He I 


2.05 


0.14 




4070 


[S II] (blend) 


1.99 


0.08 




4101 


H<5 


24.60 


0.11 




4188 


cm 


0.27 


0.05 




4200 


Hen 


0.61 


0.04 




4267 


CII 


0.78 


0.08 




4340 


H7 


43.710 


0.11 




4363 


[O III] 


10.32 


0.07 




4387 


He I 


0.39 


0.10 




4471 


He I 


3.40 


0.07 




4543 


He II 


1.22 


0.07 




4648 


cm 


2.75 


0.57 




4685.7 


He II 


41.86 


0.11 


38.1 


4711.4 


[Ar IV] + He I 


4.05 


0.06 


4.0 


4740.2 


[Ar IV] 


3.15 


0.07 


3.3 


4861.3 


H/3 


100.00 


0.00 


100.0 


4923 


He I 


1.23 


0.06 




4958.9 


[0 III] 


355.65 


0.51 


359.4 


5006.8 


[O III] 


1062.02 


189.91 


1047.2 


5412 


Hell 


3.69 


0.04 




5518 


[CI III] 


0.50 


0.04 




5539 


[CI III] 


0.46 


0.04 




5755 


[NE] 


0.30 


0.11 




5805 


CIV 


0.90 


0.11 




5876 


He I 


11.04 


0.13 




6300 


[O I] 


1.37 


0.07 




6312 


[S III] 


2.12 


0.04 




6363 


[O I] 


0.46 


0.06 




6408 


He II 


0.19 


0.04 




6548 


[Nil] 


2.01 


0.06 




6563 


Ha 


279.07 


49.76 




6584 


[Nil] 


11.53 


0.07 




6678 


He I 


2.65 


0.04 




6716 


[S II] 


2.00 


0.04 




6731 


[S II] 


2.97 


0.04 




6891 


He II 


0.35 


0.04 




7065 


He I 


2.33 


0.07 




7136 


[Ar III] 


9.36 


0.07 




7175 


He I 


0.57 


0.06 




7211 




2.98 


0.09 




7322 


[Oil] 


2.26 


0.08 




7362 




0.43 


0.07 





Observed H/3 flux = 1.29 ± 0.01 xl0~ 13 (EFOSC2) and 1.06 xlO" 13 
(ISIS) ergs cm" 2 s" 1 



comp uter code which is based on a genetic algorithm 1 Gesicki et al. 
2006). The line profiles include thermal broadening, organized (ex- 
pansion) velocities, and turbulence. However, for the modelled Sgr 
PNe we did not find evidence for a turbulent component. 

5.2 Kinematical models for three PNe 

5.2.1 He 2-436 

The spectra are spatially unresolved, as shown by the fact that the 
lines are unsplit (Fig. [3}- Four lines are suitable for the velocity 
analysis; the weakest of these ([S III] 6312A) was filtered to smooth 
the noise. The lines are approximately of Gaussian shape, with 
[O III] and [S III] slightly narrower than the other two. 

In accordance with the HST images we build an appropri- 
ate photoionization model assuming the density distribution is 
concentrated towards the inner radius. The chemical comp osition 
and line ratios hav e been adopted from Dudzia k~et alJ (2000) and 
Wals h~et"al] ll997t) . A parabola-like expansion velocity resulted in 
acceptable fit to the observed line profiles. 

The adopted solution is shown in Fig. [3] Panels on the left 
present the model structure. Shown are: the assumed radial hydro- 
gen density distribution together with the computed electron den- 
sity and temperature; the derived velocity distribution; and the ra- 
dial distributions of emissivities for lines shown in right panels, 
normalized to unity. The other panels compare models with obser- 
vations. In the middle the observed HST surface brightness profiles 
from Fig.|2|are compared to the computed ones. The observations 
(dots and dashes) are rescaled to the normalized model profiles 
(draw line). The panels on the right compare the computed (solid 
lines) with the observed (circles) emission line profiles. 

5.2.2 StWr2-21 

The 2-arcsec nebula is expected to be partially spatially resolved by 
the spectrograph slit. Only four emission lines are strong enough 
to be considered in our analysis. The He II and H/3 lines are un- 
split, [O III] shows a partial split, and [N II] is split more clearly but 
is also more asymmetric. For the model fitting we used only the 
stronger (blue) peaks from the asymmetric couple. 

The best-fit model that reproduces both the H/3 flux and the 
line ratios of Table [2] is a density-bounded one. Such a model re- 
produces the very similar size of the images in the two different fil- 
ters (Fig.0, because in this case the ionization stratification is very 
weak and the emissivity distribution of [O III] and Ha follows the 
same density distribution. (The computed emisivities show for He II 
and [N II] a clear stratification, but no images are available to test 
this.) To reproduce the small size of the central cavity, we needed 
to apply a density distribution concentrated towards the centre with 
inner radius close to the central star. 

The linearly increasing velocity field reproduces relatively 
well the four observed emission lines. 

5.2.3 Wray 16-423 

An analysis of the velocity field of t he nebula Wray 16-423 was 
published in lGesicki & Ziilstri] <2003l) . The high quality spectrum 
allowed for extraction of nine emission lines which strongly con- 
strained the deduced velocity field. The result indicated a high ve- 
locity spike located inside the nebula and also suggested a presence 
of a low density extended outer region. 

The previous modeling was done without the benefit of a 
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high-resolution image. We redid the modeling, searching for a 
density distribution that reproduces as well as possible the new 
HST imag es. The obtain ed density field is similar to that from 
iGesicki & ZiilstrJ l2003h . but the maximum is placed somewhat 
further from the centre and the outer density is higher. 

With the improved density structure we returned to the old 
spectra. The genetic algorithm found a velocity field similar to that 
from lGesicki & Ziilstral 120031) with a high velocity region inside 
the nebula. The fit co nsidered eight of the nine lines shown in 
IGesicki & ZiilstraH2003h (excluding the most noisy one) but in Fig. 
|5|we present only the most representative four lines. The velocity 
profile is more complicated than found for the other two PNe, but 
this is partly caused by the fact that more lines are available and 
that therefore the velocities as function of emissivities are better 
constrained. This allowed us to search for a higher-order solution. 
The low-order behaviour is similar to the other objects. 



Table 3. Parameters and results of the velocity field analysis 



He 2-436 Wray 16-423 StWr2-21 



log F(H/3) 


-11.4 


-11.89 


-12.6 


[O/H] 


-0.52 


-0.51 


-0.35 


Teff 


58100 


82000 


146000 


l0gT efl 


4.76 


4.92 


5.16 




3.89 


3.6 


2.95 


model R out [pc] 


0.03 


0.095 


0.16 


model ilinner IP C ] 


0.001 


0.008 


0.02 


model Mi on 


0.084 


0.40 


0.25 


outer boundary 


ioniz 


ioniz 


density 


angular diam [arcsec] 


0.5 


1.6 


2.6 


Vkver [km/s] 


19 


32 


32 


turbulence [km/s] 











dyn. age [yr] 


1650 


3100 


5200 


dyn. core mass 


0.62 


0.61 


0.61 



5.2.4 Results 

Table [3] presents the summary of the parameters of Torun models 
for the three PNe. The important parameter are the mass-weighted 
expansion velocity, Vkver, the dynamical age, and the core mass for 
which the dynamical and stellar time scales are equal. 

As expected, there is an age progression with He 2-436 be- 
ing the youngest and StWr2-21 the oldest nebula. The difference 
is smaller than the size would imply, as He 2-436 has a small ex- 
pansion velocity. The derived velocities are within the usual range 
for Galactic objects IGesicki & Ziilstra 2000). Available data sug- 
gests that th e metallicity has little effect on the PN expansion 
iGesicki & Ziilstn]l2003h . This is in contrast to AGB outfl ow ve- 
lociti es which are strongly dependent on metallicity I Marshall et al. 
2004). The reason is that the PN expansion is largely dominated by 



the overpressure within the ionized region, and has accelerated sig- 
nificantly during the post- AGB phase. The kinematic models show 
that all three stars have comparable core (final) masses. There is a 
systematic uncertainty, in that some assumptions are made on how 
the nebula accelerates, and that the radius used for the dynamical 
age is taken as 0.8 of the outer (model) radius. (The latter because 
this is approximately the mass-weighted radius.) 

All three nebulae show an asymmetry in the line profiles with 
the blue component stronger. The difference is very minor in He 2- 
436, larger in Wray 16-423 and strongest in StWr2-21. In all cases 
do the ionization levels tracing the outer-most regions (such as [O I] 
and [N II]) show the largest effect. The eff ect of the asymmetry on 
the line modelling is discussed in Gesicki & Ziilstra 1 2003): the ve- 
locity structures are more regular than the density and/or tempera- 
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ture profile, as shown by the line profiles. The cause for the blue 
asymmetry is not known, but one can speculate about a relation to 
the movement of the nebulae through the gas in the Galactic halo. 



6 DISCUSSION: STARS AND NEBULAE 
6.1 Photo-ionization modelling 

Emission lines are sensitive to the local conditions in the nebulae: 
density, temperature, radiation field, and elemental abundances. 
This sensitivity allows for accurate models, with improving ac- 
curacy when more lines are available. Deep spectroscopic obser- 
vations of Wray 16-423 and He 2-436 have previousl y been used 
to ob tain photo-ionization models of both nebulae t Dudzi ak et alJ 
2000). Neither object could be well described by a set of unique 
conditions, but each required two sectors with different covering 
factors and conditions. Their model parameters are listed in Table|4] 
The models have filling factor unity and assume constant pressure 
within each zone. Table|4]also lists the radii and thickness derived 
from the HST images. The models were derived before the HST 
images became available: the values predicted by the models are 
largely confirmed by the new images presented above. The struc- 
tures observed by HST confirm the need for at least two-zone mod- 
els. The complexity in the images exceeds what can be reproduced 
in the models, but the large-scal e geometry confi rms the predictive 
capability of the models of Dudz iak et al J 120 00). 

Photoionization modelling of StWr 2-21, again using NEBU, 
is discussed in detail in Pequignot et al. (2006, hereafter Paper II). 
Models are based on the fluxes listed in Table|2| but with minor re- 
visions suggested by a spectrum synthesis. The synthesis was used 
to investigate the contribution of weak, not individually detected, 
lines on the spectrum. It also allowed us to resolve blends, impor- 
tant at the moderate resolution of the observations, and to detect 
Wolf-Rayet features (see Section l6~3l . 

As for the two objects modelled by Dudzia k et alJ fcOOOl) . at 
least two sectors are necessary to reproduce the observed line flux 
ratios. With the constraint of the nebular size derived from the HST 
images, the two-sector models for StWr2-21 turn out to have a 
combined covering factor much less than unity. A low-density gas 
with no specific signature in the spectrum may fill in most of the 
'empty' sector, as confirmed by trial calculations. Given that the 
overwelmingly dominant sector is strongly matter-bounded, the ex- 
istence of such a component is not too surprising. Although it may 
not emit much of the line flux, this sector may encompass a very 
significant fraction of the ionized mass, which is therefore only a 
lower limit in the two-sector model. 

The question of the uniqueness of the solution is discussed in 
Paper II. It is difficult to obtain unique models in the case of StWr 2- 
21. The carbon abundance, which controls the energy output of the 
nebula and therefore the effective temperature of the central star, is 
poorly constrained for lack of UV spectra. Most optical recombi- 
nation lines of carbon have mildly accurate fluxes. Moreover, fol- 
lowing a suggestion made for other PNe, these recombination lines 
may partly trac e cool H-poor clumps embedded in the nebula (e.g. 
iLiu et all2004l:|Peauignot et all20 03). The NEBU model selected 
in Table[4]is therefore typical but not unique. Nonetheless, the oxy- 
gen abundance is relatively well defined. Abundances derived from 
this model are provided in Tableland discussed in Section l6~71 

We used a Cloudy model as an independent check on the 
NEB U results. We used version 05.10.06 of Cloudy, last described 
by iFerland et alJ 1 19981) and modified to include additional stel- 



Table 5. Absolute (dereddened) stellar magnitudes for three Sgr PNe, for 
a distance of 25 kpc. V -band values are fro m the HST data above. B-band 
values are derived from Walsh e t aljll997i) . 



PN 


E(B-V) 


V 


B 


He 2-436 


0.28 


-0.61 ±0.1 


-0.3: 


StWr 2-21 





+2.70 ±0.1 




Wray 16-423 


0.14 


+1.43 ±0.2 


+1.7: 



lar atmospheres as described below. The two codes are indepen- 
dent, and use different assumptions regarding the nebular struc- 
tures. A differen t stellar atmosphere m odel was used: whereas in 
the NEBU code a Clese & Middlemass 1 1987) model was adopted, 
in the Cloudy model we used a He/C (i.e. H-poor) model of lRauchl 
(2003). Cloudy also allows the use of dust within the ionized re- 
gion: we ran models both with and without dust. The dustless mod- 
els were marginally better, but in view of the uncertainties in the 
modelling, it seems doubtful that this is significant. Cloudy con- 
firmed the need for a two-component model to obtain a reason- 
able fit. In the Cloudy models the microscopic volume filling factor 
(Osterbrock 1989) was treated as a free parameter. It converged 
to unity, tentatively indicating that the nebular material is homoge- 
neous. But this required a much less luminous star than found in the 
NEBU model. Radii and density of the components are also some- 
what different. The oxygen lines are poorly fitted in comparison to 
the NEBU results. Convergence to a full solution with Cloudy was 
not attempted. NEBU outputs are presented in Paper II. 

6.2 The central stars: HR diagram 

The absolute stellar magnitudes, dereddened and converted to John- 
son V, are listed in Table|5| for an assumed distance of 25 kpc. The 
data are corrected for extinction, derived from the Ha/H/3 ratio. 
A standard extinction curve with Rv = 3.1 is assumed. B-band 
value are also give n: these are from the long-slit spectroscopy of 
IWalshet all i 19971) and are of lower accuracy. 

Fig. [5] shows the location of the stars on the HR diagram, 
based on the photoionization models. The uncertainty on the dis- 
tance is upto 10 per cent, leading to a systemic uncertainty in log L 
of 0.1 dex, which is not included in the error bars. The depth of the 
Sgr ma in body has not been measured but is likely less than 1 kpc 
iMaiewski et al.l2003t) . 

Th e dashed lines in Fig. |6| show the models of B lockeJ 
Jl995bl) . labelled by their core mass. The speed of evolution is a 
strong function of core mass. The 'dynamical core mass' in Table 
|3|gives the mass which is consistent with the dynamical age of the 
nebula. In all three cases, this mass is around 0.61-0.62 M©. (The 
systematic uncertainty in the distance of 10 percent, introduces a 
~ 1 per cent error in the core mass, whic h is much less tha n the 
model uncertainties of ~ 0.02 Mm : iGesicki & Ziilstral200d) . The 
photoionization models yield a stellar luminosity which is a little 
lower than corresponds to these masses. However, the models may 
underestimate the stellar luminosity if there is little ionized gas in 
some directions. 

The stellar magnitudes can also be used to derive observed 
locations in the HR diagram. We used the stellar atmosphere mod- 
els of IClegg & M iddlemass 1 1987) to calculate predicted V-band 
magnitudes as function of temperature and luminosity. The re- 
sults are shown in F ig. [6] the squares show the luminosities of 
Clees & Middlemass 1 1987) models at specific temperatures, fit- 
ted to the observed V-band magnitudes. The long-dashed lines con- 
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Table 4. Results of the multi-sector photo-ionization models. Values for StWr2-21 are from Paper II. Values for He 2-436 and Wray 16-423 are from 
Dudzi ajTet alj 1200 0') Observed radii (third and fourth rows) are from HST data presented here (except for BoBn 1). The covering factors are the fractions 
of the total solid angle, as seen from the central star, occupied by each model sector. The penultimate row gives the radial optical thickness at the ionization 
threshold of H where required. The last row gives the fraction of the total ionized mass in each sector. 





He 2-436 


Wray 16-423 


StWr2-21 


BoBn 1 


(units) 


T e ff 


0.70 10 5 


1.07 10 5 


1.19 10 5 


1.30 10 5 


K 


L 


5.40 10 3 


4.35 10 3 


4.10 10 3 


5.20 10 3 




M ion 


0.02 


0.25 


0.06 


0.27 


M© 


dinner (HST, observed) 


610" 


-3. 


210" 


-2 


410" 


-2 






pc 


AR (HST, observed) 


310 


-2 


510 


-2 


810 


-2 


(~ 


41.1) 


pc 


Model sectors 


I 


II 


I 


II 


I 


II 


I 


II 




Ionization/density bounded 


ionization 


ionization 


ionization 


density 


ionization 


density 


ionization 


density 




<N H > 


2.710 s 


2.0 10 4 


9.5 10 4 


3.6 10 4 


1.9 10 3 


1.6 10 3 


1.3 10 3 


1.1 10 3 


cm -3 


Rinn (model) 


8.1 10~ 3 


8.1 10~ 3 


2.3 10~ 2 


2.3 10" 2 


4.9 10~ 2 


4.9 10~ 2 


1.3 10~ 2 


1.3 10~ 2 


pc 


AR (model) 


7.1 10~ 4 


2.3 10" 2 


2.810" 2 


5.9 10" 2 


9.0 10~ 2 


4.1 10" 2 


1.9 10 - 1 


1.610" 1 


pc 


Covering factor 


0.62 


0.38 


0.17 


0.83 


0.0074 


0.372 


0.028 


0.437 




r(13.6eV) 








6.7 




1.2 




3.1 




Mass fraction 


0.100 


0.900 


0.108 


0.892 


0.092 


0.907 


0.115 


0.885 





nect the squares; the stars are expected to fall along these lines. 
(In fact, if the dynamical mass is assumed, the stars should be lo- 
cated on the corresponding intersection between the long-dashed 
and the solid lines.) For He 2-436, the bright stellar magnitude in- 
dicates a cooler star than given by the photoionization models. The 
most likely explanatio n is that for this H -poor [WC] star, the H+He 
models of lClegg & Middlemassl f 19871) overpredict the flux at ion- 
izing wavelengths. For Wray 16-423 a similar but smaller effect is 
seen. The stellar magnitude of StWr2-21 was used as a constraint 
in the NEBU model. The relatively bright star rules out the low lu- 
minosity of the Torun model, confirming that much of the stellar 
radiation escapes in directions with low density gas as required in 
the NEBU model. We note that the results do depen d on the atmo- 
sphere model, e.g. using the more modern models of Rauch <2003l) 
could yield different results. The composition of the stellar atmo- 
spheres, and the existence of a strong wind, introduce a significant 
additional uncertainty. 



6.3 The central stars: WR lines 

After subtracting all the nebular emission lines from the spectrum 
of StWr2-21 by means of spectrum synthesis, two broad features 
at 3820A and 5 805 A were revealed. The resulting stellar spectrum 
is discussed in Paper II. The presence of broad O VI 381 1+34A and 
C IV 5 801+ 12A emission poin ts to a very hot WR star of subclass 
[WO] lAcker & Neinerll2003l) . The intensity ratio in the spectrum 
synthesis (see paper II) is Ovi 3822/ClV 5805 « 7.8 ± 2. The 
uncertainties arise from the ill-defined 'half line width' of the 3822 
line (taken as 1300km/s compared to 750km/s for 5805), and the 
blending of the stellar 5805 line with (weaker) nebular 5805 emis- 
sion. In a difficult blend, we find C IV 4658/C IV 5805—1-1.5, and 
there is an indication for O VI 5290/C IV 5805 ~ 0.25 (where both 
lines are assumed to have the same width). These line ratios in- 
dicate a [W O 21 star, perhaps mid way between subclass 2 and 3, 
according to l Acker & Neinerl 120031) . 

The two PN e previously known also show WR lines 
I Walsh^^j^^g^We reclassified the se stars using the sc heme of 
Acker & Neiner 1 2003 f), and the data ofiwalsh et alJ <1997l) subject 
to the reanalysis of Dudziak et al. (2000). The adopted subclasses 
are listed in Table Q The main constraint is the absence of some 
important lines, suggesting they are much weaker than C IV. For 



He 2-436, the He II line points at a possible [WO 4] class. Wray 16- 
423 is more uncertain, but the detected lines are consistent with 
[WC4]. The higher temperature of the star of Wray 16-423 is not 
reflected in its subclass. However, this is not unusual among WR 
nuclei, where the subclass also traces the strength of the stellar 
wind 1 Crowther 1999). The narrow lines of Wray 16-423 could also 
be interpreted as indicative of a weak emission-line star (wels). 

All three Sgr PN with deep spectra th erefore show evidence 
for WR-type emission lines. iGornvl 1200 ll) estimates that 6.5 per 
cent of Galactic PNe show [WRl-typ e central stars. A higher frac- 
tion is found in the Galactic Bulge iGornv et alj|2004l) . predom- 
inantly with late subtypes . To this should b e added the weak 
emission-line stars (wels; iGesicki et al J 1200(11) which occur with 
roughly equal frequency as do WR-type central stars. To find all 
three studied Sgr objects to fit in this category, supports models 
where the fundamental parameters of the progenitor star (mass, 
metallicity) determine whether a star develops a WR wind. This is 
in contrast to models where the occurence is random, caused by the 
precis e timing of the last thermal pulse (for a review see lBlockeJ 
2001). There is no evidence for WR-central stars to be more mas- 
sive th an non-emission-line stars lPottaschll996l) . but lGesicki et alJ 
(2006) find that Galactic [WC] stars show a narrower, intermediate 
range of core masses, compared to non-emission-line stars. The Sgr 
PNe fit within this mass range. 



6.4 The IR-[WC] stars 

Among Galactic [WC] stars, there is a subgroup with unusually 
strong IRAS emission. He 2-436 has similar characteristics: it is as- 
sociated with an IRAS point source <Ziilstraetalll994) . and scal- 
ing its 12-micron flux (0.22 Jy) to a typical distance of galactic PNe 
of 2.5 kpc would put it among the brightest such objects. The class 
of IR-[WC] stars is discussed bv lZiilstr alf200ll) . 

The Gala ctic IR-[WC] stars have late subtypes, of 8-11 
JZiilstralEool) . He 2-436 shows a much earlier subtype. Three 
further extragalactic [WC] stars are detected by IRAS: SMP58 
and SMP61 in the Large Magellanic Cloud and SMC N6 in 
the Small Magellanic Cloud JZiilstra et alJll994tlPena et alll997t 
IStasifiska et alj 120041) (the last object has no l2-[im detection). 
The two LMC objects have subclass [WC5-6] , and SMP N6 has 
[WO 4] on the scale of Acker & Neiner] 1200 3l) . The extragalactic 
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Figure 6. Location of the three stars in the HR diagram. The solid lines are the theoretical evolutionary models o flB16ckeil fT995b). labelled with the core mass 
in solar units. The long-dashed lines show the location of the stars based on the V-band magnitude: from right to left He 2-436, Wray 16-423 and StWr2-21. 
The square symbols show the specific models of Clegg & Middlemass 1 1987) used in deriving these lines. The photoionization models are shown by the plus 
signs: in all cases they are on or to the left of the corresponding long-dashed line. 



stars have much earlier subclassses than comparable Galactic stars. 
The strong IRAS flux densities require a high dust optical depth, 
and indicate a dense, compact nebula. The cause of the differences 
between the stars in different galaxies is not clear. 



6.5 Nebular structures 

6.5.1 Bipolarity 

The HST images show a very similar global morphology between 
the three imaged objects. All show a torus-like structure, with 
lower-density regions in the polar directions. In the classification of 
lGreigHl97ll) . all three objects are of class B, indicative of enhanced 
AGB mass loss in the equatorial plane. The images do not show 
the'butterfly' morphology of the more extreme bipolar nebulae, 
suggesting the density contrast between e quator and poles is not ex- 
treme (in the scheme of Balick & Frank (2002), the objects would 
be considered as intermediate between elliptical and bipolar). The 
fact that all three nebulae show similar structures is interesting in 
view of suggestions that the primary factor determining the nebular 
morphology is stellar progenitor mass ([Corradi & Schwarz 1995). 
Revi ews on PNe morp hologies can be found in Kwok 1 2000J) and 
iBalick&Frankl 1200211 . 

The outer structures of the two larger objects are more compli- 
cated; relatively bright and showing symmetry axes which are not 
well aligned with those of the inner torus. The outer reg ions are sen- 
sitive to the mass-loss evolution on the AGB {Monrea l-Ibero et alJ 
2005) and interaction with the ISM {Warei ng et all2006l) . There is 
an indication that the asymmetry in the outer stuctures is along the 



direction of motion of the galaxy (towards WNW): interaction with 
the tenuous medium in the Galactic halo surrounding Sgr may be 
affecting the nebular expansion. 

Assuming that the higher-density tori can be associated with 
the ionization-bounded sectors in the NEBU models (Table[4), the 
tori of Wray 16-243 and He 2-436 contain roughly 10% of the total 
mass of the nebulae. This should be taken with care, as the mod- 
els are imperfect representations. For StWr2-21, the HST image 
suggests that the torus does not show a strong ionization front, con- 
sistent with the model where this sector has a very low covering 
factor. For Wray 16-423, the torus has a clear sharp edge as ex- 
pected from an ionization front. In both cases, the mass fraction 
may be underestimated. For He 2-436, the dense sector has a very 
small thickness which suggests that the ionization front has barely 
penetrated the high density inner shell. The low total mass of the 
model also suggests that a large fraction of its mass is not yet ion- 
ized. 

The masses in the nebulae are typical compared to Galactic 
objects. The bipolar structures are also typical, with half of Galac- 
tic PNe showing similar morphology (Corr adi & Schwar3ll995F) . 
He 2-436 has perhaps a more extreme structure, with a very pro- 
nounced, compact torus and a high density contrast. A number of 
Galactic post-AGB stars show similar structures, but these have 
typically cool and unevolved central stars. 

IStanghellini et alJ l2003h finds similar statistics on morphol- 
ogy between the LMC, SMC and the Galaxy. Our sample agrees 
with the prevalence of the same basic morphological types, regard- 
less of environment. This is an important result: it provides evi- 
dence against models where the (mild) bipolarity is caused by sub- 
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stellar/planetary companions iLivio & Soke3 r2002). as this would 
tend to introduce a metallicity dependence which is not observed. 
Strongly bipolar nebulae may be the product of a unique evolution 
for each object (Soker 2002): for the milder stuctures seen here, a 
single evolutionary cause may be preferred. 



6.5.2 Double shell 

The presence or absence of a double-shell structure is an impor- 
tant characteristic of PNe. The evolution of a PN is dominated by a 
few effects: (i) the expansion and density profile of the progenitor 
AGB wind; (ii) the fast wind from the central star blowing into the 
AGB wind; (iii) the ioniz ation of the swept-up n ebula and the resul- 
tant p ressure imbalance iPerin otto et aljE004t ISchonberne r"et alj 
2005). The fast wind sets up a compressed rim at the inner radius: 
this is the usual PN. The over-pressure in the ionized region sets up 
a swept-up shell further out. This leads to the 'doubl e shell' c onsist- 
ing of the bright inner rim and the fainter outer shell. lPerinotto et all 
(2004) presents and extensive set of hydrodynamical models for 
spherical PNe, to show how this structure evolves as function of 
AGB mass loss and central star mass. The double shell appears 
after the ionization front breaks out from the inner rim, and the 
nebula becomes density-bounded. In earlier phases where the PN 
is ionization bounded, only the rim is seen. 

A double shell is evident in the radial profiles of Wray 16- 
423 and StWr2-21 (Fig. |2}. The nebula can be approximated 
by the hydrodynamic al model f or a 0.605 Mq central star (se- 
quence 6 of lPerinotto et alj |2004). Here the nebula becomes den- 
sity bounded after approximately 3500 yr. The outer radius of the 
rim of Wray 16-423 (~ 1.5 x 10 17 cm), and the brightness con- 
trast between the rim and shell of a factor of 4-5, are to first or- 
der consistent with the corresponding Perinotto model in the early 
density-bounded phase. (The spherical asymmetry in the nebulae 
will cause differences in evolution with respect to the symetric 
models.) This model uses a high AGB mass-lo ss rate, peaking at 
3 x 10" 4 Mq yr" 1 (case B2 of lBlockerl fl995j) was used). It also 
has a too high progenitor mass of 3 Mq . The size of the rim in- 
creases with lower mass-loss rate, but it is not a strong function. 
The models suggests that a uniform M > 3 x 10~ 5 Mq yr -1 can 
also be accommodated in order to explain the observed structure. A 
high peak mass-loss rate is however more likely, as AGB mass-loss 
rates are known to increase with time. 

The outer shell of StWr2-21 is less distinct because of a low 
contrast: the brightness contrast is a factor of 2-3. The Perinotto 
models show an increasing contrast with time. As StWr2-21 is 
the more evolved nebula, this suggests its evolution has not been 
identical to that of Wray 16-423. However, the differences are rela- 
tively minor. The structure of the nebula changes strongly after the 
star enters the cooling track, but this has not yet happened even for 
StWr2-21 (sectionlO. 

He 2-436 shows no indication for such a double shell struc- 
ture. This is as expected in its early phase of evolution, where 
the ionization front i s still trapped inside the rim. The models of 
iDudziak et alj J2000) require a central dense inner shell together 
with a lower-density component. However, comparison with the 
images suggests this can be interpreted as the torus and the polar 
flows, rather than a radial double-shell structure. 



6.5.3 Dynamics 

The models o flPerinotto et aljfcopj) and Schonb erner et alj<2005l) 
also predict velocity fields, which can be compared to those derived 
here. The flow velocities are sensitive to the wind density structure. 

The youngest object, He 2-436, shows a decreasing velocity 
with radius over the highest density gas. As in all cases, the veloc- 
ities increase again towards the outer edge (ionization front). The 
predicted velocities in this phase show the same behaviour as ob- 
served, but the predicted velocities are around lOkrns -1 while the 
observed expansion is faster. For Wray 16-423, a complicated ve- 
locity field is found which shows only a vague resemblance to the 
predicted one (within the highest density gas, the slight decline of 
expansion velocity with radius is consistent with the models, as is 
the value of the velocities in this region of around 20kms~ . The 
observed sharp increase near the inner edge is not predicted, nor is 
the increase in the low density outer gas. 

For StWr2-21, the velocities increase linearly with radius, in 
contrast to the other two modelled nebulae. This is a predicted fea- 
ture of more evolved nebulae, and is a better indication of its evo- 
lutionary state than the density contrast of the rim and shell. The 
observed velocities (reaching 50km s -1 near the outer edge) are 
somewhat large compared to the model values. 

The issue of high observed expansion velocities is addressed 
by Schonberner et al. 1 2005), who shows that velocities are higher 
if a steeper initial density distribution, p oc r~ °, is assumed. Our 
expan sion velocities, compared to Fig. 12 of Ischonberne r et"all 
(2005), indicate values of a = 2.5 for He 2-436, and a ~ 3 for 
the other two objects. This density distribution is expected from an 
AGB wind which increases in time. 

We note that both the Torun models and the Schoenberner 
models are one-dimensional. This is a significant limitation for in- 
trinsically 3-d objects. Observed structure in the velocity field, such 
as found in Wray 16-423, may result from averaging over different 
directions. The derived parameters in the Torun models: average 
expansion velocity, dynamical age and mass, have been shown to 
be robust against assumptions regarding nebular structure, but ef- 
fect such as turbulence are not. Three-dimensional kinematic mod- 
els are being developed (e.g. iMorisset et aljEo05T) but full hydro- 
dynamical models are still restricted to 1-D. 



6.6 Evolution 

Both stars and nebulae are consistent with an evolutionary sequence 
where He-2436 is the youngest, Wray 16-423 more evolved and 
StWr2-21 is the most evolved nebula, but where the initial condi- 
tions were very similar. The dynamical ages (Table|3j clearly show 
this sequence; the nebular radius increases whilst the density drops 
and the ionization front moves out. For StWr2-21, very little of the 

mass is still neutral. 

The IRAS detection of He 2-436 jZiilstra & Walsrll996h also 
shows this object to be young with strong heating of the dust in 
proximity to the star. The absence of an IRAS detection for the 
other two objects is consistent with their much cooler expected 
dust. 

The comparison with the hydrodynamical models, especially 
based on the velocity fields, suggests that the AGB mass-loss rates 
increased with time, and reached peak values of ~ 10 -4 Mq yr -1 . 
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6.7 Abundances 

Abundances for the Sgr PNe are listed in Table |5| V alues for 
He 2-436 and Wray 16-423 are from lPudziak et al]J200fJ) . StWr2- 
21 abundances are from the models discussed in Section 16,11 
The carbon abundance, based on the C II recombination line at 
4267A, is uncertain. The Cloudy model indicated lower abun- 
dances for StWr2-21, by about 0.15 dex. We consider the fully 
converged NEBU model as more accurate, but either model shows 
that StWr2-21 has higher abundances than the other PNe. 

We also give abundances for BoBn 1. Unlike for the other Sgr 
PNe, a relatively good IUE spectrum is available for this PN. Abun- 
dances of BoBn 1 have been obtained by H oward et alUl997h from 
models. However, for consistency with the other three PNe, the 
abundances derived from a two-sector NEBU model worked out 
by Pequignot (unpublished) are quoted in Table |6| The N, O and 
Ne abunda nces are in excelle nt agreement with a recent compila- 
tion by Perinotto et al ] j2004h (who do not provide S and Ar). The 
NEBU results on the other hand differ somewhat (but for different 
reasons) from either Howard et al. ( 1997) or the re-compilation by 
the same group iHenrv et al. 2004). 

Table |6] sho ws the eleme ntal abundances relative to the solar 
abundance o f lLoddersl f2003\ These abundances have been under 
discussion, and a downward revision to [O/H] = 8.66 has been pro- 
posed l Asplund et al. 2004). This revision has proble ms reproduc- 
ing th e known depth of the solar convection zone iBasu & Antial 

The two previously analyzed PNe have identical abundances, 
which agree well with those of the carbon stars of Sgr. StWr2- 
21 has a higher abundance. BoBn 1 has a complicated abundance 
pattern. This is an extraordinary object in that the neon abundance 
exceeds that of oxygen. Such an abundance pattern can only be 
caused by dredge-up of nuclear-burning products. The neon en- 
hancement may indicate that the star has suffered a Very late 
Ther mal Pulse (VL TP), similar to for example Sakurai's Object 
lHaiduk et all2005h . Alternatively, binary coalescence can be con- 
sidered. C,N,0 and Ne cannot be used to constrain its progenitor 
abundance: their high abundances reflect their formation within the 
star (Peauignot et al. 2000; Peauignot & Tsamis 2005). However, 
both the sulfur and argon abundances of the PN indicate that [Fe/H] 
< — 2 and possibly as low as —2.3, considering that, for low metal - 
licit y stars, [Fe/Hl tends to be a few 0. 1 dex smaller than [Ar/H] 
(e.g. lCavrel et all2004h . 

The [S/O] of the three PNe in the core of Sgr is -1.60 for 
StWr2-21, and —1.70 and —1.77 for the other two objects. This is 
close to but perhaps a litt le lower than the solar value listed in Table 
|6|(— 1.59: the revi sion oflA splund et al. 1 2004) would give an even 
higher solar ratio). Caffau et al. 1 2005) present sulphur abundances 
for five stars in Terzan 7, a globular cluster in Sgr with a metallicity 
similar to He 2-436 and Wray 16-423. They find [S/O] = -1.71, 
identical to our result. Sulphur and oxygen are both produced in 
massive supernovae, but the lower ratio in the lower metallicity ob- 
jects may be due to primary oxygen production in low-metallicity 
AGB stars l Peaui gnot et all2 000). 



7 DISCUSSION: THE SGR GALAXY 
7.1 Sgr Luminosity 

We can estimate the luminosity of the Sgr galaxy from its PNe pop- 
ulation. The specific fre quency of PNe is fa irly well determined 
within the Local Group: Magrini "e7aT] J2003I) find log iV(PN) = 



log Lv — 6.92, in solar units. This predicts for Sgr a value of 
Lv = 10 7 ' 5 L0, or Afv = —14.1, provided the PNe population in 
Sgr is now reasonably complete. The lack of PNe near the centre of 
Sgr may raise do ubt on this, however th e centre is also lacking in 
globular clusters. Maiewski et al. (2003) derive Afv = —13.27 for 
the main body, assuming a King profile. We confirm their suspicion 
that this underestimates the total luminosity. 

Assuming a ma ss-to-light ratio of 20-50 llbat a et ail 1 19971 
iMaiewski et alj l2004). the total mass of Sgr becomes M = 1— 
1.5 x 10 9 Mq. This is close to the minimum value required for 
its survival up to the present jlbata et alj[l997l) . The fact that one 
PN is found in the tail makes it likely that the tails account for a 
significant fraction of the Sgr mass, but not a dominant fraction. 

7.2 The enrichment history of Sgr 

StWr2-21 shows noticeably high abundances compared to He 2- 
436 and Wray 16-423. The identical abundances of the latter sug- 
gest their progenitor stars formed in a single, major star-formation 
burst which took place in a well-mixed ISM. The burst has been 
dated at 5 Gyr ago, based on the AGB colour, and breaking the 
age-m etallicity degeneracy with the PN abundances ( Dudziak et al. 
2000). 

Evidence for a m ore metal-rich population was presented by 
Bonif acio et alJ<2004l> . based on spectroscopy of 10 red-giant stars. 
They derive [Fe/H]~ —0.25. This is in agreement with StWr2-21, 
for which we find [O/H] = —0.23. The high metallicity suggests a 
recent star formation event. The last major bu rst of star f ormation in 
Sgr o c cured between . 5 and 3 Gyr ago (e.g. lLavden & Saraiedinil 
200Q): lBonifacio et alj J2004l) date it at 1 Gyr or less. (More than 
one event may have occured.) Assuming 1 and 5 Gyr, progenitor 
masses for the PNe would be 2.2 and 1.3Mq, respectively. 

The fact that only one of the PNe shows this high m etallicity 
does not agree with the finding of iBonifacio et alj J2004) that the 
large majority of their stars shows the higher metallicity. Allowance 
should be made for the small number of objects in both surveys. 
The difference could also be related to the location of their stars: 
IBonifacio et"ai]l2004h surveyed a region located midway between 
the core region of Sgr and the three PNe. Bell azzini et alj 1 2006) 
find that the large majority of the Sgr population has a metallicity 
in the range [M/H] = -0.4 to -0.7 and an age of 8.0 ± 1.5 Gyr. 
The more metal-rich, younger population traced by StWr2-21 may 
not be uniformly distributed. 

All four PNe are strongly carbon rich. The morphology 
of the Sgr dwarf ha s been traced e xtensively with M giants 
iMaiewski et alj|2003t iBellazzini et alj|2006l) . However, the PNe 
confirm that even for the most metal-rich and youngest population, 
3rd dredge-up causes the AGB stars to evolve into carbon stars. 

BoBn 1 derives from a much more metal-poor population. The 
existence and significance of such a population in Sgr is shown by 
the metal-poor clusters M54 (which m ay be the nucleus of Sgr) an d 
Arp2 have metallicity [Fe/H] = -1.8 <Lavden & SaraiedinilpOOOf) . 
and Ter 8 is reported as [Fe/H] = —2.0 I Lavden & Saraiedini 1997). 
The RR Lyrae population also shows a min or but significant popu- 
lation with [Fe/H]Ss —2 iCseresnies 2001). BoBn 1, if its associa- 
tion with Sgr is confirmed, may derive from this earliest population 
of Sgr. 

7.3 Evolution of Sgr 

The abundance distributions of Dwarf Spheroidal galaxies, in- 
cluding Sgr, differ from dwarf irregular galaxies of similar lumi- 
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Table 6. Comparison of abundances, given on a logarithmic scale where H = 12. 



Elem. 


Wray 16-423 


He 2-436 


StWr2-21 


BoBnl 


Galactic 
PNe a 


Solar 6 © 


Wray - © 


He 2-© 


St-0 


BB-© 


GPNe - © 


H 


12.00±.00 


12.00 


12.00 


12.00±.00 


12.00±.00 


12.00 












He 


11.03±.01 


11.03±.01 


11.00 


11.01 


11.05±.03 


10.99 


+0.04 


+0.04 


+0.01 


+0.02 


+0.06 


C 


8.86±.06 


9.06±.09 


9.00: 


8.83 


8.81±.30 


8.46 


+0.40 


+0.60 


+0.54 


+0.37 


+0.35 


N 


7.68±.05 


7.42±.06 


7.88 


7.84 


8.14±.20 


7.90 


-0.22 


-0.48 


-0.02 


-0.06 


+0.24 


O 


8.33±.02 


8.36±.06 


8.53 


7.75 


8.69±.15 


8.76 


-0.43 


-0.40 


-0.23 


-1.01 


-0.07 


Ne 


7.55±.03 


7.54±.06 


7.72 


7.83 


8.10±.15 


7.95 


-0.40 


-0.41 


-0.23 


-0.12 


+0.15 


Mg 


6.98±.30 










7.62 


-0.64: 










S 


6.67±.04 


6.59±.05 


6.93 


5.16 


6.91±.30 


7.26 


-0.59 


-0.67 


-0.33 


-2.10 


-0.35 


CI 


4.89±.18 










5.33 


-0.44: 










Ar 


5.95±.07 


5.78±.08 


6.10 


4.46 


6.38±.30 


6.62 


-0.67 


-0.84 


-0.52 


-2.16 


-0.24 


K 


4.65±.22 










5.18 


-0.53: 











a Mean composition and scatter for non-Type I Gal actic PNe iKinssbursh & Barlow 1994) 
Solar system abundances are from lLodderl<2003l) 



nosities. Both contain old stars of very low abundances (which 
in dwarf irregulars s how a more spheroidal-like distribution: 
iMinniti & Ziilstrall99rj) . But whereas younger stars in dlrr galax- 
ies still are at low abun dances, with a st rong relation between [O/H] 
and galaxy luminosity (Pilvugin 2001), for dSph galaxies younger 
stars are more metal-rich I Grebel et al. 2003). For galaxies of simi- 
lar mass and luminosities, the dSph galaxies have a dominant pop- 
ulatio n whic h is much more metal-rich than do dlrr galaxies (e.g. 
lMateoll998l) . Different reasons have been suggested: more efficient 
enrichment, less infall of unprocessed gas or a large amount of tidal 
stripping I Grebel et al. 2003). dlrr galaxies may have low efficiency 
in enrichment due to blow-out of stellar ejecta by supernovae, 
or they may benefit from continuous infall of unprocessed extra- 
galactic HI. dSph galaxies may have lacked any subsequent gas in- 
fall (Lavden & Saraiedini 1 2000) argue for a closed-box chemical 
evolution in Sgr). Present dSph systems may also have had a much 
higher luminosity/mass in the past with the present galaxy being 
only a remnant of the original system. It is however likely that the 
cause is related to the location of the galaxies: dwarf spheroidals are 
found near large galaxies while dwarf irregulars are more isolated. 

The ISM in which the recent star formation of Sgr occured 
is unlikely to have been original. The orbit of Sgr takes it through 
the Galactic plane, outside the solar circle, every lGyr or so. This 
is believed to be a possible cause of sudden star formation bursts, 
but the passage will also strip Sgr of its gas. At present there is no 
evidence for ISM within the main body of Sgr although stripped HI 
has been found along its orbit iPutman et al. 2004). 

Two sources of ISM renewal can be considered. The first is 
mass loss from stars, as traced by the planetary nebulae. The three 
PNe in the main body of Sgr suggest a current stellar death rate of 
about 10~ 3 yr _1 . For progenitor masses of 1.3-1.5 Mq, and stel- 
lar remnants of 0.6 Mq, the mass return rate to the ISM becomes 
5 x 10 -4 Mq yr _1 . Over 1 Gyr, the approximate orbital period of 
Sgr, the ISM can reach ~ 10 6 Mq which may be sufficient to trig- 
ger star formation. In the absence of significant rotation in the stel- 
lar population, this lost mass will congregate in the potential well. 
Thus, the next star formation eve nt will occur near the centre of the 
galaxy, as is observed in Fornax ISaviane et all20 00). As the orig- 
inal low metallicity ISM has been preferentially removed, the en- 
richment efficiency will be considerably higher than in the closed- 



box model. All these points are consistent with the properties of the 
dSph galaxies, including Sgr. 

A second possiblity to be considered is that the dSph captured 
gas from the Galactic halo. Assuming that the halo contains a reser- 
voir of cooling gas originating from a Galactic Fountain, a galaxy 
in orbit within the halo can be expected to capture some of this 
gas. This can explain the very high metallicities, as the halo gas has 
metallicities appropriate for the host galaxy. In combination with 
the self-enrichment scenario above, this would give an even faster 
increase of metallicity with time. It is noticeable that the metallicity 
of the most recent star formation event in Sgr has a metallicity, as 
shown by StWr2-21, not below that of regions of the outer Galactic 
plane. 

We finally note the similarity with the Fornax dwarf, a galaxy 
with similar mass and properties to Sgr. It also shows evidence 
for an extended star formation history with increasing metallic- 
ity. Both the old and the youngest population (at ~ 2 10 8 yr even 
younger than Sgr) have similar metal licity to the ir counterparts in 
Sgr, [Fe/H] = -0.4, to -2.0 iSaviane et all200d) . Fornax has suf- 
fered much less ram pressure stripping and tidal stripping from the 
Galaxy, being in a much more distant orbit. This suggests that the 
fast chemical evolution is not dominated by gas removal through 
ram pressure. Instead gas accretion, or lack thereof, may drive the 
evolution. Enrichment in Fornax and Sgr may have been acceler- 
ated by a lack of an external source of metal-poor gas (as compared 
to isolated dlrr galaxies), or by accretion of enriched gas expelled 
by our Galaxy. 

The origin of the difference in enrichment history between 
dSph and dlrr galaxies is still an enigma ) Greb eTet al I2003T) . It is 
worth testing whether the isolated transition galaxies between dlrr 
and dSph, such as Phoenix, happen to be located in regions of a low 
density of accretable, non-enriched gas. 



7.4 The tidal tails 

The tidal tails of Sgr encircle the Milky Way Galaxy in a compli- 
cated pattern. The s ystem clearly has had a long history of disrup- 
tion. The models of lLaw et alj y005) show how different parts of 
the tails trace different distortion events. 

The Sgr tidal tails are an important, but not d ominant, c ontri- 
bution to the Galactic halo. On statistical grounds. iBellazzini et alj 
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1 2003) argue that 4 of the 35 halo globular clusters should be linked 
with Sgr debris. Since approximately 10 halo PNe are known, our 
finding that one of these (BoBn 1) is a tidal member of Sgr is in line 
with the globula r cluster fra ction. 

iBellazzini et alJ 1200 3l) do not consider the globular clusters 
beyond 60kpc, as falling beyond the apogee of the Sgr orbit. But 
the tidal streams are found out to 80-100 kpc lLaw et al]|2005h 
and at least one of their discarded distant clusters (NGC2419) is 
attributed to Sgr debris ( Newbe rg et al J 120031) . No halo PNe are 
known at anywhere near such distances, possibly because of faint- 
ness. 

The tidal tails trace the older , metal-poor population of Sgr 
(Alard 2001 ; Bellazzini et al. 2002), as shown by its RR Lyrae stars 
(Cseresnies 2001). IBellazzini et al. (2006) argue that the dominant 
metal-rich population should also be distributed around the tidal 
tails.) The globular clusters indicate that the [Fe/H] of the tidal tails 
is around —1.5 to —1.9. This population becomes dominant in the 
tails because of the age-metallicity gradient: the earlier disruptions 
affected the older, less metal-rich population. 

Within the Local Group, populations with [Fe/H] ^ — 1 show 
indications for a deficit of PNe IZiilstr afcOollMagrini et all2003l) . 
related to a lower mass-loss efficiency on the AGB. This would 
enhance the fraction of peculiar objects (e.g. VLTP nebulae, for 
which the formation rate does not depend on metallicity), to which 
category BoBn 1 likely belongs. 



8 CONCLUSIONS 

The main points of our study of the PNe population of the Sagittar- 
ius dwarf spheroidal galaxy can be summarized as follows. 

1. Two new PNe are discovered in Sgr, bringing the total num- 
ber to four. One is located in the leading tidal tail, the other three 
are located in the main body, south of the main centre of Sgr. The 
PNe population (assuming the sample is complete) indicates an ab- 
solute visual magnitude of My — —14.1 for Sgr, including its 
leading and trailing tails. A stellar death rate of 10 -3 per year is 
derived. 

2. HST images of the three PNe in the main body of Sgr show 
very similar, bipolar morphologies, consisting of a denser equato- 
rial torus and lower density, more extended polar regions. This type 
of morphology is also common among Galactic PNe, and therefore 
the origin is unlikely to be dependent on metallicity-dependent ef- 
fects, as would be expected if planetary companions provide the an- 
gular momentum required for the equatorially-enhanced mass loss. 
Whereas for strongly bipolar PNe it has been argued that every ob- 
ject is the product of a unique evolution <Soker 2002). for the mildly 
bipolar nebulae a single cause is preferred. 

3. Analysis of the central stars shows that all are on the post- 
AGB track, prior to entering the cooling track. The effective tem- 
peratures are between 7x 10 4 and 1.1 x 10 5 K, and masses are close 
to 0.61 M . Post-AGB ages are between 1000 and 5000 years. 

4. All three analyzed PNe show WR emission lines, with early 
subclasses of [WO 2] to [WC4], WR-type central stars of PNe 
form a minority population among Galactic PNe. Their prominence 
in the Sgr galaxy suggests that the progenitor star characteristics 
(mass, metallicity) determine the occurence of such stars. This is 
in contrast to models where the (random) timing of the last thermal 
pulse is the crucial factor. 

5. Photoionization models for the three PNe with HST images 
yield abundances of [O/H] = —0.23 for one object and [O/H] = 
—0.55 for the other two. The first object confirms the reality of 



the m etal-rich population uncovered recently by iBonifacio et all 
(2004), and suggests it accounts for a substantial part of the Sgr 
stellar population, albeit not its majority. The identical abundances 
of the other two objects argue for a short-lived burst in a well-mixed 
ISM: this burst has been dated at 5 Gyr ago. The higher metallicity 
object may date from a burst roughly 1 Gyr ago. The fourth object, 
BoBn 1, located in the tidal tail, has very unusual abundances and a 
non-standard history, but a progenitor metallicity of [Fe/H]~ —2 is 
indicated. This object may belong to the oldest population of Sgr. 

6. The mass-weighted expansion velocities are similar to those 
of Galactic PNe, and show no indication of a metallicity effect. 
Kinematic ages and core masses are derived: all three objects have 
similar core masses of around 0.61 Mq. The age-velocity relation, 
compared to hydrodynamical models, suggest a steep initial density 
distribution, of p oc r -3 . This indicates that the preceding AGB 
stars showed mass-loss rates increasing with time. 

7. The ISM reformation rate of Sgr due to stellar mass loss is 
5 x 10~ 4 Mq yr^ 1 . This amounts to 10 6 Mq over one orbit, which 
may be sufficient to trigger renewed star formation. 

8. The steep metallicity-age gradient, and the very high re- 
cent abundance, must be due to a combination of stripping of older 
ISM from Sgr, with enrichment from its own stellar mass loss. The 
fact that Fornax, which is in a much wider orbit around the Milky 
Way, shows very similar recent enrichment, may indicate that ram- 
pressure stripping is not the dominant effect. We suggest that accre- 
tion of gas originating from the Milky Way may play a role. This is 
based on the fact that the most metal-rich population has a metal- 
licity not too different from that of the outer regions of the Milky 
Way Galaxy. 

9. We finally note that He 2-436 provides the sole direct detec- 
tion of dust in a dwarf spheroidal galaxy, to date. 
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